Amicable pairs and aliquot cycles for elliptic curves 
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Abstract. An amicable pair for an elliptic curve E/Q is a pair 
of primes {p,q) of good reduction for E satisfying ^Ep{¥p) — q 
and ^Eq{¥q) = p. In this paper we study elliptic amicable pairs 
and analogously defined longer elliptic aliquot cycles. We show 
that there exist elliptic curves with arbitrarily long aliqout cycles, 
but that CM elliptic curves (with j ^ 0) have no aliqout cycles 
of length greater than two. We give conjectural formulas for the 
frequency of amicable pairs. For CM curves, the derivation of pre- 
cise conjectural formulas involves a detailed analysis of the values 
of the Crossencharacter evaluated at primes p in End(i?) having 
the property that #£'p(Fp) is prime. This is especially intricate 
for the family of curves with J = 0. 



Introduction 

Let -E/Q be an elliptic curve. In this paper we study pairs of 
primes (p, g) such that E has good reduction at p and q and such 
that the reductions -Ep and Eq of E at p and q satisfy 

#Ep(Fp) = g and #E,(Fg) = p. 

By analogy with a classical problem in number theory (cf. Remark [7]), 
we call {p,q) an amicable pair for the elliptic curve E/Q. 

Example 1. Searching for amicable pairs using primes smaller than 10^ 
on the two elliptic curves 

El : + y = — X and E2 : y'^ + y = x^ + x^, 

yields one amicable pair on the curve Ei, 

(1622311,1622471), 
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and four amicable pairs on the curve E2, 

(853,883), (77761,77999), (1147339,1148359), (1447429,1447561). 
Example 2. The curve 

: = + 2 

exhibits strikingly different amicable pair behavior. There are more 
than 800 amicable pairs for using primes smaller that 10^, the first 
few of which are 

(13, 19), (139, 163), (541, 571), (613, 661), (757, 787), (1693, 1741). 

One objective of this note is to present theoretical and numerical 
evidence for the following conjecture. 

Conjecture 3. Let E/Q be an elliptic curve, let 

Qe {X) = # { amicable pairs {p, q) for E/Q with p < q and p < 

be the amicable pair counting function, and assume that there are in- 
finitely many primes p such that ^Ep{¥p) is prime. 

(a) If E does not have complex multiplication, then 



QsiX) X — as X ^00, 

{\ogXy 

where the implied constants depend on E. 
(b) If E has complex multiplication, then there is a constant Ae > 
such that 

Qe{X) ~ Aej-^. 

[logxy 

We do not believe that it is clear, a priori, why there should be such 
a striking difference between the CM and the non-CM cases. We first 
discovered this phenomenon experimentally; subsequently we found an 
explanation based on Theorem [131 which says that if E/Q has CM and 
if g = 4I^Ep{¥p) is prime, then there are generally only two possible 
values for ^Eg{¥q), one of which is p. (The situation for j{E) = is 
considerably more complicated; see Section [61 ) This contrasts with the 
non-CM case, where ^Eg(¥g) seems to be free to range throughout the 
Hasse interval. We refer the reader to Conjectures [TSl and [101 for more 
precise versions of the CM part of Conjecture [3l 

The frequency of primes p such that jj^Ep{¥p) is prime or almost 
prime has been studied by a number of authors. In Section[Tlwe discuss 
what is known and what is conjectured concerning this problem. 
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Generalizing the notion of amicable pair, we define an aliquot cycle 
of length £ for i^/Q to be a sequence of distinct primes {pi,p2, ■ ■ ■ ,Pe) 
such that E has good reduction at every pi and such that 

Example 4. The elliptic curve y"^ = — 25x — 8 has the aliquot triple 
(83, 79, 73). The elliptic curve 

E:y'^ = x^ + 176209333661915432764478X + 60625229794681596832262 

has an ahquot cycle (23, 31, 41, 47, 59, 67, 73, 79, 71, 61, 53, 43, 37, 29) of 
length 14. 

In Section [3] we give an heuristic argument suggesting that the count- 
ing function for aliquot cycles of length i for non-CM elliptic curves 
grows like V^/(logX)^. The rough idea is to assume that if g = 
^Ep(¥p) is prime, then the trace values ag{E) = q + 1 — j^Eq(¥g) are 
(more-or-less) equidistributed within the appropriate Hasse interval. 

In Section m we give an elementary construction (Theorem [TTj) using 
the prime number theorem, the Chinese remainder theorem, and a 
result of Deuring, to prove that for every i there exists an elliptic 
curve E/Q with an aliquot cycle of length i. 

We next consider the case of elliptic curves having complex multi- 
plication. These curves exhibit strikingly different behavior from their 
non-CM counterparts. Our first result (Theorem [T3l) says that if E/Q 
has CM with j{E) ^ 0, and if g = ^Ep{¥p) is prime, then there are 
only two possible values for ^Eg{¥q), namely p and 2q + 2—p. Assum- 
ing each is equally likely (which seems to be the case experimentally), 
this explains why CM curves have so many amicable pairs. The proof 
involves first proving that p and q split in End(£^), and then relating 
the values of the Grossencharacter at primes lying above p and q. 

Theorem [T3] can also be used to show that a CM curve with j 7^ 
has no aliquot cycles of length 3 or greater; see Corollary [TH This 
stands in contrast to Theorem [TTl which says that there exist curves 
with arbitrarily long aliquot cycles. 

We finally turn to the j = curves y"^ = x^ + k, whose complicated 
analysis is given in a lengthy Section [61 For prime values of k, we give a 
precise conjectural formula for the counting function of amicable pairs 
that depends on the value of k modulo 36. For example, if k is prime 
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and /c = 1 or 19 (mod 36), then we conjecture that 

#|p < X : p is part of an amicable pair) 1 1 
lim — ^ 7 = ^ = - + , (1) 

x-*^ < X : #^p(Fp) is prime} 6 3A; - 9 

while if A; = 11 or 23 (mod 36), then the limiting value in ([1]) is (con- 
jecturally) equal to | + -^^zq- There are similar formulas for the other 
congruence classes. 

The derivation of these formulas is in two parts. First, by analyz- 
ing the values of the Grossencharacter and using sextic reciprocity, we 
prove that (p, q) is an amicable pair if and only if (^-^^) g ^^^) g = 1- 
If the values of ipsip) modulo k were equidistributed as p varies, we 
would conjecture that the number of amicable pairs is governed by the 
proportion of A G O/kO satisfying (^^x^)6 ~ (Here O = 'End(E) = 
Z[(l + y/—3)/2].) This is almost true, but the allowable values of A are 
often restricted by further conditions on {j)^- Sorting out these restric- 
tions gives a precise conjectural value for the limit ([1]) in terms of the 
sizes of certain subsets of O/kO. 

The second part of the proof is to derive explicit formulas for the 
sizes of these sets. This is done by relating the points in these sets to 
the (9/fc(9-points on a certain family of curves C*-'^'''^ of genus four. We 
count these points by explicitly decomposing the Jacobian of C'^^'^'^ into 
a product of four j = elliptic curves and using the Grossencharacter 
formula for the number of points on such curves. The resulting for- 
mulas are quite involved, especially in the case that k splits in O, but 
eventually most of the terms cancel, leaving a relatively compact for- 
mula. We have no good explanation for why the final formula has such 
a simple form; see Remark [38] for a discussion of the delicacy of the 
computation. 

The conjectures in this paper are supported by heuristic arguments 
and, especially for CM curves, by theorems describing the allowable 
values of the Grossencharacter ipE- But heuristic arguments have been 
known to fail, and indeed our CM argument depends on the assumption 
that V'-b(P) mod k is uniformly distributed among its allowable values, 
where we claim to have characterized the set of allowable values. It 
is thus reassuring that extensive experiments are in close agreement 
with the conjectural values derived by theory. These experiments are 
described in Section [71 

Finally, in Section [8] we explain where we first ran across amica- 
ble pairs and aliquot cycles for elliptic curves, and we describe some 
possible generalizations that deserve further study. 
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1. How OFTEN IS i^Ep(¥p) PRIME? 

If an elliptic curve E/Q is to have any amicable pairs or aliquot 
cycles, then it is clearly necessary that there exist primes p such that 
Ep(¥p) is prime. The question of existence and density of such primes 
has been studied by various authors. 

Remark 5. If i?(Q)tors 7^ {O}, then ^Ep(¥p) will be composite for all 
but finitely many p, since -E(Q)tors ^ Epi¥p) for all p \ 2Ae/q- Using 
this observation, it is quite easy to produce curves having no nontrivial 
aliquot cycles, for example, the curves y"^ = + x and = + 1. 

More generally, there may be a local obstruction associated with 
the representation Gal(Q/Q) Aut(i?tors) that forces ^Ep(¥p) to be 
composite for all but finitely many p; see [22]. We quote a special case 
of a conjecture of Koblitz, as modified by Zywina. 

Conjecture 6. (Koblitz [9], Zywina [29]) Let E/Q be an elliptic curve, 
and let 

Me^X) = ^ {^primes p < X such that ^Ep(¥p) is prime^ 

count how often E modulo p has a prime number of points. Then there 
is a constant Ce/q such that 

Me{X) ~ Ce,^^^,. 

Further, Ce/q > if and only if there are infinitely many primes p 
such that ^Ep(¥p) is prime. 

Koblitz and Zywina give formulas for the constant Ce/q in terms 
of the image of the representation Gal(Q/Q) — >■ Aut(£'tors)- In prin- 
ciple this allows one to approximate Ce/q to high precision, and they 
give a number of examples. For additional material on the probability 
that ^Ep(¥p) is prime or almost prime, see [2] and [S]. 

2. Aliquot cycles and amicable pairs for elliptic curves 

We formally give the following definitions as previously described in 
the introduction. 

Definition. Let E'/Q be an elliptic curve. An aliquot cycle of length i 
for i?/Q is a sequence of distinct primes {pi,p2, . . . ,pe) such that E has 
good reduction at every pi and such that 

#Ep,_A^p,_,) = Pi, #Ep^{¥pJ=p,. 
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An aliquot cycle is normalized if pi = minpj. Every aliquot cycle can 
be normalized by a cyclic shift of its elements. An amicable pair is an 
aliquot cycle of length two. 

Remark 7. Classically, an amicable pair is a pair of integers (m, n) 
satisfying a{m) = n and &{n) = m, where 5"(?2) = cr(n) — n is the 
sum of the proper divisors of n. Similarly, a number n is perfect if 
5"(n) = n, and a (classical) aliquot cycle is a list of distinct integers 
(ni,n2, . . .,ni) satisfying 

a{ni) = n2, a{n2) = n^, ... a{ni-i) = n^, a{ni) = ni. 

(Numbers appearing in an aliquot cycle are also called sociable num- 
bers.) Perfect numbers and amicable pairs were studied in ancient 
Greece, and aliquot cycles of all lengths continue to attract interest to 
the present day. See for example [H [5], [251 l2H] ■ 

By analogy with the classical case, one might call an aliquot cycle (p) 
of length one for E/Q a perfect prime, but such primes have already 
been given a name. They are called anomalous primes and appear 
as exceptional cases in diverse applications; see for example [121 fT3] . 
In particular, anomalous primes are to be avoided in cryptography 
because the elliptic curve discrete logarithm problem (ECDLP) for 
anomalous primes can be solved in linear time fi6[ \T7\ [22] . 

We begin our study of aliquot cycles with the following general ob- 
servation concerning amicable pairs. 

Proposition 8. Let E/Q be an elliptic curve, and let {p,q) be a nor- 
malized amicable pair for E/Q with p > 5. Then 



Proof. The fact that p is odd and q = ^Ep(¥p) = p + 1 — ap is prime 
implies in particular that Op 7^ 0, so has ordinary reduction at p. 
(This is where we use the assumption that p > 5; cf. [201 Exercise 5.10].) 
Reversing the roles of p and q shows that E also has ordinary reduction 
at q. 

The assumption that {p, q) is an amicable pair is equivalent to the 
assertions that 



End(Ep/Fp) ® Q = End(Eg/Fg) ® Q. 



q = p+l- ap 



and 



p = q+l- a, 



and then a little bit of algebra shows that 



Op - 4p = - 4g. 



(2) 
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The field End(i?p/Fp) ® Q is generated by the Frobenius element 
Frobp(a;) = x^, which is a root of 

-apT + p = 0. 

Thus 

End(Ep/Fp) ® Q ^ Q (y/a2 - Ap ) . 

The analogous formula is true for q, and then ([2]) completes the proof 
of the proposition. □ 

3. Counting aliquot cycles for non-CM elliptic curves 

In this section we study the ahquot cycle counting function 

n ( Y\ — 41- f normalized aliquot cycles {pi, ■ ■ ■ ,pe) 
JeA^) - # length i for E/Q satisfying pi < X 

Conjecture 9. Let E/Q be an elliptic curve that does not have com- 
plex multiplication, and assume that there are infinitely many primes p 
such that ^Ep(¥p) is prime. Then the aliquot cycle counting function 
satisfies 



y X 
(logX)* 

where the implied positive constants depend on E and i, but are inde- 
pendent of X . 

Remark 10. As noted in Section [HI an aliquot cycle (p) of length one 
consists of a single anomalous prime. In this case, Conjecture [9] follows 
from a general conjecture of Lang and Trotter [IQl, which predicts the 
stronger result Qs,i(X) ~ c^/X/\ogX. 

We give an heuristic argument in support of Conjecture M To ease 
notation, let 

Np = #4(Fp). 
Then, setting pi = p, we have 

Prob(p is part of an aliquot cycle of length i) 



Prob 



~ ( n Prob(pi+i =^ Np,^ is prime) j Prob(A/p^ =Pi)- (3) 

(We ignore the small probabilty that there is some i < i such that iVp. 
is equal to an earlier pj.) 



P2 == Np-^^ is prime and =^ Np^ is prime and 
. . . and Pi =^ Xp^_^ is prime and Np^ = pi 
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Under our assumption that Np is prime for infinitely many p, Con- 
jecture E] says that 

Prob(A^„ is prime) ■; , 

logp 

and since 

Pi+i = Np^ =Pi + 0{^), 

every term in the sequence p = pi,p2, ■ ■ ■ ,pe satisfies Pi = p + 0{y/p). 
Hence 

ProbfAL is prime) ~ . 

log Pi logp 

Substituting this into ([3]) gives 

-n 1 /^p is part of an ahquot\ 1 n \ 

[ cycle of length £ J " " ^^^^^^^ = P'^' 

In order to estimate the last factor, we use the Sato-Tate conjec- 
ture [20l C.21.1], which says that as q varies, the values of Ng are 
distributed in the interval [q + 1 — 2y/q, q + 1 + 1^fq\ according to the 
Sato-Tate distribution, 

#|g<X:a< ^^^~^'^ <&l^ 7r(X) ■ - [' VT~^ dt. 

[ 2^ J TT ^ 

(See |24j for a proof of the Sato-Tate conjecture in certain cases, al- 
though our use of the conjecture is purely heuristic.) Then for primes p 
and q = p + 0{y/p), a rough estimate gives 



PToh{Ng = p) X — -. (5) 



1 

Combining (jlj) and ([5]) yields 



Prob P^^^ aliquot\ 1 

\^ cycle of length £ ) v^(logp)^"^' 



We now estimate the number of normalized aliquot cycles of length 
whose initial prime is less than X. 

O iX) ~ Prob initial element of a normalized\ 

eA J ~ 2-^ y aliquot cycle of length £ J 



X V r^, 

^v^(logp)^ 
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It only remains to use the rough approximation 
^/(X)^ J2 f{n\ogn)^ / f{t\ogt)dt^ / f{u)f^ 

p<X n<X/logX J OgU 

to obtain 

^ .... I du VX 



V^(logM)^ ^ logM (logX)^ 



4. Aliquot cycles of arbitrary length 

Theorem 11. For every i > 1 there exists an elliptic curve E/Q 
that has an aliquot cycle of length i. More generally, for any positive 
integers ii, . . . ,ir there exists an elliptic curve E/Q that has distinct 
aliquot cycles of lengths £i, . . . ,ir- 

Proof. A theorem of Deuring [3] (vastly generalized by Waterhouse [27], 
see also Riick [15]) says that if p is a prime and t is an integer satisfy- 
ing |t| < 2y/p, then there exists an elliptic curve E/¥p satisfying 

#E(Fp)=p + l-t 

In other words, every Frobenius trace in the Hasse interval for p actually 
occurs as the trace of an elliptic curve defined over ¥p. 

Now fix i and let pi,p2, . . . ,pe be a sequence of primes with the 
property that 

\Pi + l- Pi+i\ < 2v^ for all 1 < 2 < £, (6) 

where by convention we set p^+i = pi- It is easy enough to find such a 
sequence. To be precise, we can use a weak form of the prime number 
theorem pH, Theorem 4.7] that says that there are positive constants a 
and b such that the n^^ prime g„ satisfies 

an log(n) < gri < 6nlog(r;,). 

It follows that for any given i, if we choose n to be sufficiently large, 
then 

This implies that the sequence of primes g„+2, • • • , '?n+^) satis- 

fies (E]), so we take this to be our sequence {pi, . . . ,pi). 

Applying the theorem of Deuring cited earlier, for each pi we can 
find an elliptic curve Ei/¥p- satisfying 
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(This includes the case i = £, in which case pe+i = pi-) We now use 
the Chinese remainder theorem on the coefficients of the Weierstrass 
equations for Ei, Ei to find an elhptic curve E/Q satisfying 

E mod Pi = Ei for all 1 < i < £. 

Then by construction, the sequence (pi, . . . ,pe) is an ahquot cycle of 
length i for E/Q. 

In a similar fashion, we can construct elliptic curves over Q that have 
aliquot cycles of any specified lengths using different sets of primes, and 
then we can Chinese remainder these curves to obtain a single elliptic 
curve over Q with any specified number of aliquot cycles of any specified 
lengths. □ 

Remark 12. The algorithm described in Theorem [TT] works well in prac- 
tice, although it naturally yields equations having large coefficients. 
We used it in Example S] to find an aliquot cycle of length 14. Here's 
another example. The following elliptic curve has an aliquot cycle of 
length 25, starting with the prime p = 41. 

2^2 = 3,3 _,_ 4545482133607498579268567738514832922289740324532X 

+ 595867265462112118291430245894379464967885794713. 



5. Amicable pairs for CM curves with j 7^ 

Our next goal is to formulate and provide evidence for more precise 
versions of the CM part of Conjecture O A key observation is that 
if E has CM, then the assumption that q = jj^Ep{¥p) is prime severely 
limits the possible values of Eq{¥q). It turns out that the case of elliptic 
curves with j{E) = is significantly more complicated than the other 
cases, so we deal with the j{E) 7^ curves in this section and leave the 
j{E) = curves for the next section. 

Theorem 13. Let E/Q be an elliptic curve and assume: 
(1) E has complex multiplication by an order O in a quadratic imag- 



inary field K = Q{V-D) . 
(2) p and q are primes of good reduction for E with p > 5 and 

q = i^EpiWpj 



l+^ 



(3) j{E) 7^ 0, or equivalently, (9 7^ Z 
Then D = 3 (mod 4), and either 

#4(F,)=p or #4(F,) = 2g + 2-p. 
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Theorem [13] has an interesting consequence concerning the allowable 
lengths of aliquot cycles for CM elliptic curves. This may be compared 
with Theorem HI which says that there exist (necessarily non-CM) 
curves having aliqout cycles of arbitrary length, and with Conjecture [H 
which implies that every non-CM elliptic curve has aliqout cycles of ar- 
bitrary length provided that there are infinitely many primes p such 
that ^Ep(¥p) is prime. 

Corollary 14. A CM elliptic curve E/Q with j{E) ^ has no aliquot 
cycles of length i > 3 consisting of primes p > 5. 

Remark 15. There are various ways in which one might generalize The- 
orem [131 For example, replacing assumption (2) by the assumption 
that L is an integer such that the quantity 

g = L2- (p + l-#4(Fp))L+p 

is prime and splits in K leads to the following conclusion: 

a,iE) = ±{apiE) + 2L). 

Theorem [13] is the case L = 1. We omit the proof of the generalization, 
since it is similar and is not required in this paper. 

Remark 16. Corollary [TJ] omits curves with j{E) = 0. It turns out 
that j = curves possess a rich and complicated amicable pair struc- 
ture which will be investigated in detail in Section [6] Corollary [2l] 
gives an analogue of Theorem [13] saying that there are (often) six pos- 
sible values for ^Eq{¥q), rather than only the two possibilities given 
in Theorem [121 Using this result, we are able to prove by a detailed 
case-by-case analysis that j = curves cannot have aliquot cycles of 
length three; see Appendix [X] But we do not have a proof that there 
are no aliquot cycles of length greater than three. 

Before commencing the proofs of Theorem [121 and Corollary [TH we 
prove a basic result concerning the splitting of primes in CM fields. 

Lemma 17. Let E/Q be an elliptic curve with complex multiplication 
by K , let p > 5 be a prime of good reduction for E/Q, and suppose 
that ^Ep{¥p) is odd. Then p splits in K . 

Proof. We have 

#Ep{¥p) =p+l-ap, 

so the assumptions that p 7^ 2 and ^Ep(¥p) is odd imply that Op is odd, 
so in particular Op 7^ 0. Hence E has ordinary reduction at p. (Note 
that our assumption that p > 5 and Basse's bound \ap\ < 2y^ imply 
that p \ Qpii and only if = 0.) It follows that the field K is isomorphic 
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to End(-E'p) (S> Q, which is generated by a root of the characteristic 
polynomial — a^T + p of Frobenius. Therefore K = Q{^yap — Ap ) , 
and 

P = 

either sphts or is ramified in K. But we can rule out the latter case by 
noting that 

p ramified in K =^ p\ ^P =^ P I o-p =^ <^p = 0. 
This contradicts the fact that is odd, and hence p splits in K. □ 

Proof of Theorem [23 Up to Q- isomorphism, there are 13 elliptic curves 
defined over Q that have complex multiplication. For a list, see for 
example [T^, A §3]. There are three isomorphism classes whose con- 
ductor A'"^; is a power of two: 

E : = + X, Ne = 2^ 

E : = x^ + 6x^ + x, Ne = 2^ 

E -.y^ = x^ + 4x^ + 2x, Ne = 2l 

All three of these curves have a nontrivial two-torsion point, as do all of 
their Q/Q twists, so ^E{¥p) is even for all p > 3. Hence none of these 
curves admit an amicable pair; cf. Remark[5l The remaining CM curves 
have complex multiplication by a field Q(a/— -D) with D = 3 (mod 4). 



The endomorphism ring of E is an order in the field K = Q(v — -D) 
where D = 3 (mod 4), so it has the form 

End(^) ^ C = Z + /Z 

for some integer / > 1, which is called the conductor of O. In particu- 
lar, we have O* = {±1}, since our assumption that j{E) ^ excludes 
the case (D, /) = (3, 1). 

The theory of complex multiplication says that there is a Grossen- 
character ipE such that for every prime ideal p of Ok of residue char- 
acteristic p > 5 at which E has good reduction, we have 

(i) ^e{^) G O with i)E{^)OK = p. 

(ii) #4(Fp) = Nk/q{p) + 1 - Tt{MP))- 
See, for example, [HI Proposition 4.1] or |T9l II §9]. (Note that our 
assumption that p has residue characteristic p > 5 implies that p does 
not divide the conductor of O, since our assmption that O has class 
number one implies that the conductor of O is at most 3.) 
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We are given that p > 5 and that jj^Ep{¥p) = g is prime. It follows 
from Lemma [T7] that p splits in K, say 



pOk = pp. 



Then Fp = Fp, so 



q = #4(F,) = #4(Fp) = N^/q(1 - ^^(p)). (7) 

Notice that this implies, in particular, that q splits in K. So writing 
qOx = qq, we have 

q = NK/Q{M^))- (8) 
Comparing ([7]) and ([HD, and using the fact that i'sip) and V'-E(q) are 
in O, we see that there is a unit u & O* such that either 



M^) = u{^-Mp)) or M(])=u{1-Mp))- (9) 

(This follows from the fact that the factorization of the ideal qO is 
unique, up to switching the factors.) 

As noted earlier, we have O* = {±1}, so 

Tr(^i5(q)) = ±Tr(l - ^^(p)) from © with u = ±1, 

= ±(2-Tr(^s(p))) linearity, 

= ±(2 - (p + 1 - g)) since #4(Fp) = q, 

= ±(g + l-p). 

Hence 

#E,(F,) = #^,(Fq) = g + 1 - Tt{M^)) =q+l±iq + l-p). 
This completes the proof of Theorem [131 D 

Proof of Corollary^T^ Let (pi,p2, • • • ,Pe) be an aliquot cycle of length 
£ > 3 for -E/Q such that all Pi > 5. Since the primes in the cycle are 
distinct, Theorem [T3] tells us that 



Pi = 2pi_i + 2 - pi^2 hi 3 < i < i. 

Further, since the term in the aliquot sequence following pe is pi, we 
have 

Pi = 2p, + 2-p,_,. (10) 

Consider the linear recursion 

Ai = p, A2 = g, Ai = 2Ai_i + 2 - Ai„2 for i > 3. 

A simple calcuation shows that the general term of this recursion is 
given by the formula 

Ai = {i- l)g -{i- 2)p 2). (11) 
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Hence the right-hand side of ( fTOi) . which corresponds to pi+i, is equal 
to 

Equating this with pi, rearranging terms, and dividing by i, yields 

Pl = P2 + ^ - I- 

The same argument applied to the aliquot cycle 

iPi,Pi+u ■ ■ ■ ,Pe,PuP2, ■ ■ • 

obtained by cyclically permuting the terms in the original cycle yields 

Pi = Pi+i + i — 1 for all 1 < i < £, 

where we set p^+i = pi. Since i > 1, this shows that pi > p^+i (strict 
inequality). Hence 

Pi> P2> P3> ■■■ >Pi> Pe+1 = Pi- 
This contradiction completes the proof of Corollary [TH □ 

We now use Theorem [13] to give an heuristic justification for the 
following conjecture. 

Conjecture 18. Let E/Q be an elliptic curve with with complex mul- 
tiplication, and assume that j{E) ^ 0. Define counting functions 

Ae(X) = {^primes p < X such that ^Ep(¥p) is prime^, 

Qe{X) = ^ {^amicable pairs {p,q) for E/Q with p < q and p < X}. 
Then either Me {X) is bounded, or else 

lim^ = l 

X^ocAfEiX) 4' 

We note that Conjecture [6] says that if Ne{,X) is unbounded, then it 
is asymptotic to Ce/qX/ (logX)^. So the combination of Conjectures [6] 
and [18] gives a strengthened version of the CM part of Conjecture [S] 

Our justification for Conjecture [18] is to observe that Theorem [T3] 
says that if ^Ep(¥p) = q is prime, then there are two possibilities for 
j^Eq{¥g), one of which is p. Experiments indicate that each possibility 
occurs with equal probability, and we have no theoretical reasons for 
expecting otherwise, so we will accept the hypothesis that 

Prob(#^5(F,) = p I #4(Fp) = g is prime) = ^. 

Further, if we assume Conjecture [6] then 

^fEix) 



Fioh(^^Ep(¥p) is prime | p < X) 



niX) 
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Combining these estimates yields 
#{p < X : #4(Fp) = g is prime and #^q(F,) = p] 

?a ^2 Pi'ob(#Eg(Fg) = p and ^Ep{¥p) = q is prime) 

p<X 

^ J2 Prob(#Eg(F,) = p I #4(Fp) = g is prime) 
^-^ X Prob(#£'p(Fp) is prime) 

_ Me{X) 
2 

Finally, we need to divide by 2, because Qe{X) only counts amicable 
pairs (p, g) that are normalized to satisfy p < q. 

6. Amicable pairs for CM curves with j = 

In this section we study elliptic curves having j-invariant zero. The 
analysis of amicable pairs on these curves is significantly more compli- 
cated than on all other CM elliptic curves, due primarily to the extra 
units in the endomorphism ring. In particular, experiments described 
in Section [7| suggest that the limiting value of Qe{X) /Ne{X) for the 
curve y"^ = x'^ + k varies for different values of k] see Conjecture HOI 

We continue with the Grossencharacter notation from the previous 
section and set some additional notation that will remain in effect for 
this section. We let 

^ = ^"^^ ^ ir = Q(v^), OK = nu:i 

so is a primitive sixth root of unity. We note that the unit group 
{Ok I^Ok)* is a group of order 6, and that the natural map 

^,^ = 0*^^ {Ok/^OkY 

is an isomorphism. Further, for any prime ideal p of Ok that is rela- 
tively prime to 3 and any a e Ok \ p, we recall that the sextic residue 
symbol is defined by the conditions 

G/Xe and = a^(NK/oP-i) (j^od p). 

Theorem 19. Let k E^L he a nonzero integer, let E/Q be the elliptic 
curve 

E : y'^ = + k, 
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SO E has CM by Ok, and let ipE be the Grossencharacter associated 
to E. Suppose that p > 5 and q > 5 are primes of good reduction for E 
such that 

i^EM = q. 

(a) The prime p splits in K, say pOx = PP, and satisfies 

'iIje{p){1-MP)) = 1 (mod 30k). 

(b) The ideal defined by 

q = (l - i!E{p))OK satisfies qOx = qq. 

In particular, the prime q splits in K . 

(c) The values of the Grossencharacter at p and q are related by 

1-MP)={^ i^A V^Mq). (12) 



(d) Let e G {±1}- Then the trace ag{E) = q + 1 — ^Eg{¥q) satisfies 

Remark 20. The expressions in (c) and (d) appear naturally in the 
course of proving Theorem [T9|, but we note that they may be simplified 
using Proposition [23| which says that = 1- This allows us to 

rewrite (fT2!) and (fT3|) as 

i-Mp) = (^^^^(^^^Mq), m) 

a,{E) = ±{q + 1 - p) ^ G)^©^^^^- ™ 

Proof. The fact that p splits in Ok follows from Lemma [TTl which 
proves the first part of (a). Next, as noted during the proof of Theo- 
rem [T31 the Grossencharacter of a CM elliptic curve satisfies 

#4(Fp) = ^k/q{MP)) + 1 - TrK/Q{MP))- 

Using the given value q = ^Ep(¥p), this can be written as 

q = ^K/Q{l-i^E{p))- 

Hence q = (l — i{je{c[))Ok satisfies qOx = qq, which proves (b). 

Further, both ipsip) and 1 — iPe{p) have norms that are relatively 
prime to 3. This implies first that iPe{P) = (mod 3) for some j G Z, 
and second that j is odd, since otherwise 1 — would be divisible 
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by -\/— 3. On the other hand, for any odd value of j it is easy to check 
that 

(1 -cj^>^' = 1 (mod 30k), 

so we find that 

MP) (1 - MP)) = 1 (mod 30k). (14) 

This proves the second assertion in (a). 

For the proof of (c) we use the exphcit formula for the Grossencharac- 
ter of curves of the form y"^ = + k in terms of sextic residue symbols. 
This formula says that iPe{p) = ~(^)g ^) where the generator vr is a 
primary generator for p, i.e., vr = 2 (mod 30k)- (See [3, Chapter 18, 
Theorem 4 and Section 7] or Proposition 4.1].) Reducing this 
formula for ipE modulo 3 and applying it to both of the primes p and q, 
we obtain 

f4:k\ ~\ , . , , . . /4r 



Mp)=yY) (i^odSO^) and ijE{^)=y—j {mod 30k). 

' (15) 

By definition, the ideal q is generated by 1 — iPe{p)- On the other 
hand, the Grossencharacter has the property that V's(q) generates 
the ideal q. It follows that there is a unit u G = /^g such that 
1 - iPe{P) = M^£;(q)- Using ([H]) and ([15]), we find that 

1 - MP) ^ 1 ^ (^\ (^\ 
t^Eiq) ^E{p)^E{q) VPyeVq/e ° 

Since a sixth root of unity is determined by its residue modulo 3, this 
last congruence is an equality, which completes the proof of (c). 

Using the defining property of the Grossencharacter and formula (fT2l) 
from (c), we have 



ttq 



Similarly, using the assumption that ^Ep{¥p) = q, we find that 
Trs/Q(l - MP)) = 2 - TiE/QiMP)) = 2 - (p + 1 - g) = g + 1 - p. 
Hence for e G {±1}, we have 

ag{E) = e{q + 1 - p) ^ 

Trx/QU(^y)^ (y)^ i^-MP))] =^^e/q{1-MP))- 
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We now use the following lemma, which may be applied because the 
quantity N£;/q(1 —iPe{p)) = q is neither a square nor 3 times a square. 
The lemma allows us to conclude that 

which completes the proof of (e). □ 

Lemma 21. Let a G Ok have the property the ^k/q{(^) is neither a 
square nor 3 times a square. Then 

Trx/Q(C«) = Tri^/Q(a) with ( e fi^ <^=^ C = 1- 

Proof. We have 

Trx/Q(Ca) = Trx/Q(a) Tr;^/Q((C - l)a) = 

<^=^ (C — l)a = c\^^ for some c G Z, 



<^=^ C = I oi a = c ^ 

C-1 

(Note that c is in Z because ( and a are in Ok = Z[ti;].) Suppose 
that ( ^ 1. We observe that as ( ranges over Hq \ {1}, the quan- 
tity -y/— 3/ {( — 1) takes on the five values 

1 

— UJ,1 — LU, — — CO, —UJ, —1 — UJ 

The norms of these five numbers form the set {1, 3, |}, so the norm of a 
would have the form c^, 3c^, or 3(c/2)^, contradicting the assumption 
on N^-zofa). □ 



We can use Theorem [T^ to show that for some curves with j{E) = 0, 
the conclusion of Theorem [13] is true, i.e., there are only two possible 
values for ^Eq(¥g). 

Corollary 22. Let (i G Z fee a nonzero integer, and let E be the elliptic 
curve E : y'^ = + 2d^. Let p be a prime with p \ 6d such that 
q = ^Ep{¥p) is also prime and satisfies q \ 6d. Then 

#4(F,)=p or #4(F,) = 2g + 2-p. 

Proof. Using notation from Theorem [T9l we have k = 2d^, so 

(^±\ = f2^y' = ±1 and m = = ±1. 

It follows from Theorem [TW d) that aq{E) = ±(g + 1 — p). □ 
We next prove two useful facts. 
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Proposition 23. Let k, E, p, q, and q he as in the statement of 
Theorem UR 

la) 1 — 1 = cu or u!^ 



P 

(/n (b), (— )q denotes the usual quadratic residue symbol in Z.) 

Proof, (a) If is a square modulo p, then Ep{¥p) has a nontrivial 
3-torsion point, so ^Ep{¥p) cannot be prime. Similarly, if is a 
cube modulo p, then Ep(¥p) has a nontrivial 2-torsion point, so again 
^Ep(¥p) cannot prime. Hence 

This means that (^) cannot equal 1, tu^, cu^, or to^, so it must be 
either u or u^. 

(b) We first note that for any a,/? G Ok with gcd(6,/?) = 1, we have 

5)r^(i)>(i):(i)>(i).(i), 

If, in addition, a G Z, then |Tl Chapter 18, Section 7, Lemma 2] says 
that (f) = (^^) . 

In order to prove (b), we use cubic reciprocity [7, Chapter 9, Sec- 
tion 3]. We recall that an element a G Ok is said to be primary if 
a = 2 (mod 30k)- Since iPe{P) is relatively prime to 3, there is a 
(unique) sixth root of unity C ^ Me such that (iPe{P) is primary. It 
follows from Theorem [TW a) that C~^(l ~ '^-b(P)) is also primary, and 
of course, the number 2 is primary. Hence cubic reciprocity yields 

2\ /2 



2 



CMP)J3\C-'i^-MP)) 
CMp)\ fC-\i-MP)) 



^1:e{p){i~Mp)) 



(17) 
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The primes i'Eip) and 1 — 'ipEip) are relatively prime to 2, so i'Eip) 
is congruent to either uj or 1 + u modulo 2. (Note that OkI'^Ok = 
{0, 1, cij, 1 + cu}.) Hence 

i^Eip) (1 - i^Eip)) = + o;) = 1 (mod 20k). 

Substituting into <^ shows that (J) 3(^)3 = 1- Using ([HD and its 
accompanying remark, we find that 

2\"V2\"^ f2\ f2\ f2\ f2\ f2\ (2 



vP/6 Vq/6 \P) 2\^) 2\P) 2\V ^ KpJqKQ, ^ 
which completes the proof of (b). □ 

Corollary 24. Let E/Q, p, and q be as in the statement of Theo- 
remlTR 

(a) There exists an integer A satisfying 

^2 ^ 2pg + 2p + 2q-p^ - -1 ^^^^ 
3 

(b) The trace aq{E) = g + 1 — jj^Eq{¥q) equals one of the following six 
values: 

±{q + l-p)±3A 
±{q + l-p). ^ . (19) 

Remark 25. The six possible values of jj^Eq{¥q) described in Corol- 
lary [24l^b) are ^Eq'^\¥q) for the sextic twists of Eq corresponding to 
the elements of i/^ (Gal(F ^/F^), Aut(^g)) ^ H^{Ga\(¥q/¥q), n^) = 

f;/(f;)^ 

Remark 26. Using Corollary [Ml and by-case analysis, we prove in 

Appendix |X] that j = elliptic curves have no aliquot cycles of length 
three. 

Proof, (a) We know that Tr('0£;(p)) = ap{E), so writing V'-b(P) as an 
element of Ok = ^[t^], it has the form 

,Ije{p) = +-4v^ g^^^ ^ ^ ^20) 

Since we also know that ^k/q{'^e{P)) = P, "we find that 

M^ni^^,, (21) 

Finally, the assumption that jj^Ep{¥p) = g is equivalent to ap{E) = 
p + l — q. Substituting this value into (I2T1) . a little bit of algebra shows 
that A has the form specified by (fTSl) . 
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(b) Applying ( |T2l) from Theorem [191 we find that 

TTK/Q{Mn)) = Tr/,/Q(C(l - MP))) 

for some C G Me- Using the value of iPe{P) from fl20|) with the substi- 
tution ap{E) = p + 1 — q yields 

' q + 1 — p — A\f^ 



Trx/Q(V^i?(q)) = Tr^/Q (^C 

Substituting in each of the six possible values ( & Hq and taking the 
trace yields the six values listed in flTUl) . □ 

Definition. Fix a non E : = + k. We let Mk denote the set 

{primes p > 5 of good reduction for E 
such that q = ^Ep{¥p) is also 
a prime of good reduction for E 

(This differs slightly from our earlier notation in that we are now ex- 
cluding a few primes, but this does not affect our asymptotic formulas.) 
We define a subset of Mk by 

AT™ = {p G A4 : a,{E) = ±(g + 1 - p)}, 

and we say that the primes in Xl^^ are of Type 1 for E. We write A4(X) 
for the number of primes in A4 that are less than X, and similarly 
forAAW(X). 

Only Type 1 primes can be amicable, and based on experiments, we 
expect that about half of the Type 1 primes will be part of an amicable 
pair. Let 

Qk{X) = ij^[p < X : p < q and (p, q) is an amicable pair for E^, 

i-e., Qk{X) is the number of normalized amicable pairs (p, g) on E 
with p < X. Then we have the following conjecture, where the con- 
jectured limit is |, rather than |, because Qk{X) counts amicable 

pairs {p,q) with p < q, while J\f^\x) counts both {p,q) and 

Conjecture 27. With notation as above, the proportion of Type 1 
primes that are part of a normalized amicable pair is given by 

lim — f-j = -. 

Thus in order to understand the distribution of amicable pairs on E, 
we need to study the density of the Type 1 primes in A^. 
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Remark 28. According to Corollary [2H there are six possible values 
for aq{E), two of which give Type 1 primes, so one might expect A/"]^' 
to have density | inside A/1-. This turns out not to be the case. At the 

extreme end, Corollary [22] says that A/'g^i = A/'2rf3 for any nonzero d G Z. 
The rest of this section is devoted to developing tools for calculating a 
conjectural value for limx^oo A/"]^^ (X) / J\fk{X) . This value depends on k 
in quite a complicated way; see Conjecture [301 For precise formulas 
when k is prime, see Conjecture [iQl which says that the limit should 
equal | + R{k), where R{k) is a rational function of k that depends 
on k modulo 36. 



Definition. We set the notation 

pr 

n = a (mod m) <^=^ p = a (mod m) for every prime p \ n. 
Further, for any ideal ^ C Ok we define 

Ok 



O 



AG-^:gcd(A(l-A),.^)=l 



If = kOK is principal, we write simply Of^ j^. 

Now let k E Z satisfy gcd(6, k) = 1. We define a set Mk that depends 
on k modulo 4 and on the primes dividing k modulo 9. 



pv 

(a) /c = 1 (mod 4) and k = ±1 (mod 9) 



M, 




-1 and 



(b) A; = 1 (mod 4) and k ^ ±1 (mod 9) 



M, 



A e Ol, : 



pr 

(c) = 3 (mod 4) and k = ±1 (mod 9) 



M, 



\eO 



K,k 



(d) = 3 (mod 4) and A; p ±1 (mod 9) 



M, = Ol,. 

Further, for every k we define a subset of Mk by 



M, 



[1] 



A G M. 



A(l-A) 
k 



1 
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k 


#ol, 

II K,k 


#A4 


" fe 




37 (a.l) 


1225 


408 


144 


^ = 0.3529 

17 


17 (a.2) 


287 


96 


36 


1 = 0.3750 


13 (b.l) 


121 


60 


20 


i = 0.3333 


5 (b.2) 


23 


12 


4 


i = 0.3333 


19 (c.l) 


289 


192 


72 


1 = 0.3750 


71 (c.2) 


5039 


3360 


1152 


if = 0.3429 


7(d.l) 


25 


25 


13 


if = 0.5200 


11 (d.2) 


119 


119 


47 


^ = 0.3950 



Table 1. The sets O^^^, Mk and Mj^ 



Remark 29. It is easy to check that k = ±1 (mod 9) if and only if every 
cube root of unity in Ox/kOK is itself a cube. For example, suppose 
that /c G Z is prime. If A; = — 1 (mod 9), then Ox/kOx — is a 
finite field with fc^ elements, and fig C ¥^2. Similarly, if = 1 (mod 9), 
then Ox/kOK — Ffc x IF^., and /Xg C F^. Thus in both cases, every 
cube root of unity in Ox/kOK is itself a cube. 

Conjecture 30. Let k E 2, be an integer satisfying gcd(6, /c) = 1. 
Then 

N,(X) §M, ■ '''' 

Remark 31. For small values of k it is not difficult to compute the 
sets A^fc and J^^l\ thereby obtaining an explicit (conjectural) value 
for the limit fl22|) . Table [1] gives some examples corresponding to the 
four cases (a)-(d) used to define M.k, further divided according to the 
value of k modulo 3. (The notation {x.n) after each value of k indicates 
the case x = (a),. . . ,(d) and the congruence class k = n (mod 3).) 

Our justification for Conjecture [30] uses the following weak form of 
quadratic and cubic reciprocity for the field Q(a;). 

Lemma 32. Let k E Z satisfy gcd(/c, 6) = 1, and let X G Ij[uj] satisfy 
gcd(6fc,A) = 1. 

(a) (Quadratic Reciprocity in Q(c<j)) 
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(b) (Cubic Reciprocity in Q(to')) Let ( & be the unique cube root of 
unity such that 



Then 



(X = ±1 (mod 30k) 



a), \k)Xk 



Proof. Let G Z[uj] satisfy gcd(a,/9) = gcd(a/3, 6) = 1. We start 
with the sextic reciprocity law for Q{uj) as stated in [HI Theorem 7.10]. 
This says that if a and /3 are "i^^-primary" (see [TT] for terminology), 
then 



Let p = uj"^ denote a primitive cube root of unity. Then for a G Ok 
satisfying gcd(6, a) = 1, we have by definition 



a is ii^-primary 



a = ±1 (mod 3) and 

= A + Bp with A + B = 1 (mod 4). 



(This is a corrected version of [IT], Lemma 7.9], which omits the a = 
±1 (mod 3) condition and includes a superfluous 3 | B requirement.) 
We note that if a = ±1 (mod 3), then exactly one of ±a is £^-primary. 

We now consider k E Z and A G Ok as in the statement of the 
lemma. Since k is an integer and satisfies gcd(6, k) = 1, we have 

k is E'-primary <^=^ k = ±1 (mod 3) and k^ = 1 (mod 4) 
<^=^ k = 1 (mod 4), 

so (— l)'^^^^)/^fc is ii^-primary. We also note that for any a G Ok satis- 
fying gcd(6, a) = 1, Euler's formula says that 

Zl] =(-1)^ (modaO^), (24) 

and since both sides of (IMl) are sixth roots of unity, the congruence (12^ 
is an equality. In particular, 

— 1\ / , N(fc)-1 ^ , 

-j^ = (-l)— = (-!)— = 1. (25) 

It is an easy exercise to verify that there is a unique C G /J.3 such 
that CA = ±1 (mod 3), cf. [7, Chapter 9, Proposition 9.3.5]. Then one 
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of ±(^A is i?-primary, so we can apply ( 123|) to the i?-primary numbers 
a = {-lY'^-^y^k and p = ±CA. Then 1^ becomes 



2 2 



(The second equahty comes from the fact that k"^ = 1 (mod 4).) Hence 



Using ([21]) and gives 



A;\ , , N(A)-i fc-i / / 

^ ' 1)^-^ i: r • (26) 



A/g \kj Q\k 

(We note in particular that the sign used to ensure that ^A is i?-primary 
turns out to be irrelevant because (-p)g = 1-) Cubing (126!) and us- 
ing C'^ = 1 gives the quadratic reciprocity formula in (a), and similarly 
squaring fl26l) gives the cubic reciprocity formula in (b). □ 

Justification for Conjecture Let p G A4, so Theorem [TWa) tells 
us that p splits in Ok, say pOx = PP- As in that theorem, we let 
q = (l — 'iPe{p))Ok- Then squaring Theorem [TWd) yields 



Further, Proposition [23] implies that (|) (-) =1, so we find that 



3 



The (prime) ideals p and q are generated, respectively, by the ele- 
ments ipE^p) and 1—iPe{P), and TheoremfTWa) says that these elements 
satisfy 

MP) (1 - Mp)) = 1 (mod 30k). (27) 
Hence if we choose ^ G /Xg to satisfy 

CMp) = ±1 (mod 30k), 
then (127]) says that we also have 

r'{l-MP))=±'^ (mod 30k). 

This allow us to apply cubic reciprocity (Lemma [32r b). or [71 Chapter 9, 
Section 3, Theorem 1]) to compute 

/k\ ( k \ ( k 



P) 3 \V 3 \iMP)OKj 3 V^Hl - MP))Ok 
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Hence 



We now consider how the values V's(p) are distributed in OkI^Ok 
as p varies in A4. If V were chosen completely randomly, subject only 
to p = 1 (mod 3), then we might expect the values of '^e{P) to be uni- 
formly distributed among the congruence classes in OkI^Ok- How- 
ever, Proposition [211(a) tells us that (^)^ equals either uj or cu^, i.e., it 
is a primitive sixth root of unity. Equivalently, 

and Q)^^^=o..', ,29) 

i.e., neither (^)^ nor (^)^ equals 1. This gives a constraint on the values 

of V'£;(p) for V ^ -Mc- Discarding finitely many elements of A/^, we may 
assume that p f 6A;, and then reciprocity (Lemma [32]) tells us that 

where C ^ Ms satisfies CV'b(P) = ±1 (mod 3). (Note that N(?/'£;(p)) = 
p.) Hence the constraints (129|) on '?/'_b(p) from Proposition [2^a) become 



/u / 2 \^ / 3 \ ^ 

(30) 

We now make the following two assumptions, which are supported 
by experiments: 

• For p G A/fc, the value of p mod 4 is equally likely to be 1 or 3. 

• For p G A4, the value of C in fl30l) is equally likely to be 1, cu^. 



4 

or cij . 



These assumptions have the following consequences: 

• If A; = 3 (mod 4), then the first equation in fl30l) has no effect 
on the value of V's(p) mod k. 

• If /c ^ ±1 (mod 9), i.e., if cube roots of unity in Ox/kOx are 
not necessarily cubes, then the second equation in (130|) has no 
effect on the value of iPe{P) mod k. 
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On the other hand, if A; = 1 (mod 4), then the first equation in (l30l) gives 
the constraint (^-7:^)2 = ^^^^1 similarly, if /c = ±1 (mod 9), then 
the second equation in fl30l) imposes the condition ( ^^^ '^'' ) 3 7^ 1- Thus 
considering the four cases, we see that iPe{P) is in the set M^. Further, 
we note that fl28l) says that p G A^^^ if and only if iPe{p) ^ -^fc'- Hence 
it is reasonable to conjecture that the density of AT™ in A4 is given by 
the ratio #A<W/#-^fe- □ 

Conjecture[30]is reasonably satisfactory in that the sets and A^j.^' 
are easy to compute for any particular (not-too-large) value of k. In 
the remainder of this section we derive explicit formulas for i^Aik 
and #A4fc' when k is prime. We do this by breaking them up into 
subsets of the following sort. For any ideal R C Ok and any roots of 
unity C ^ Me ^ ^ Ms? we define 



A«C) = {a.01,,,:(A) =c} 



6 

MN(C,0 = {AeA«0:(^i^)^ = {}. 

As before, if ^ = kOx is principal, we write Mk{C) and mI^\(,C,). 
Further, if S" C /Xg is any set of roots of unity, we write Mj^^S) for 
the union of Msi{() with ( ^ S. With this notation, the four cases 
defining are given by 

(a) Mk = Mk{{uj,u^}), (b) Mk = Mk{{uj,uj^uj^), 

(c) Mk = Mk{{u,u^u\uj'>}), (d) Mk = MM, 

and in all cases, m|.^^ = Aij^^S,!), where C /Xg is the set for the 
appropriate case. 

We now restrict attention to the case that G Z is a rational prime 
with gcd(6, /c) = 1. If = 2 (mod 3), so k is inert in K, then the 
computation of and m],^' takes place in the field Ox/kOx — '^k^ 
with k"^ elements. On the other hand, if /c = 1 (mod 3), so A; splits as 
kOK = then 

Ok Ok Ok 



kOK ^ ^ 
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In this case a condition such as 7^ 1 becomes more comphcated, 
since there are many ways for the product (^3(^3 to be different 
from 1. 



Proposition 33. Let k > 5 be a rational prime. The following table 
gives the values of ^Mk{S) for various subsets S (Z Hq, divided into 
cases according to whether k is split or inert in K — Q(-\/— 3 ). 





k = 1 (mod 3) 


k = 2 (mod 3) 


(a) 




l{k-l){k-3) 




(b) 




i(A;-l)(fc-3) 


\{k--l) 


(c) 


#Mfc({a;,a;^o^^o;5}) 


|(A;-l)(A;-3) 




(d) 




{k - 2f 


^-2 



Proof. Suppose first that k is inert, so Ox/kOx = ^k'^- Then ^Mi^{^q) 
simply counts the A G ¥1^2 such that A and 1 — A arc units, so is 
equal to k"^ — 2. Next, #Mfe({a;, a;^, a;^}) counts the quadratic non- 
residues in Fjfc2, of which there are \{k'^ ~ 1)- (Here the condition that 
1 — A be a unit is irrelevant, since 1 is a quadratic residue). Similarly, 
^Mk{^{uj,uj'^ ,uj^ ,uj^y) counts cubic non-residues in Ffc2, of which there 
are |(A;^ — 1). Finally, #Mfc({u;, o;^}) counts the elements that are 
neither quadratic nor cubic residues, of which there arc ^{k"^ — 1). 

Next suppose that k splits, so Ox/kOK = x F^. Then ij^Mk{^i^) 
counts (a, h) G F^ with neither a nor h equal to or 1. This gives k — 2 
possibilities for each of a and 6, so 7^Mjt(/Xg) = {k — 2)^. 

The required calculations for each of the remaining cases have much 
in common, so we will only illustrate the case of Mfc({a;,a;'^,ci;^}). Ex- 
actly I of invertible (a, b) are quadratic non-residues. Therefore, there 
are \{k—lY such elements. Of these, there are \{k— 1) of the form (1, b) 
and \{k — 1) of the form (a, 1). The set Mfc({a;, cu^, o;^}) counts invert- 
ible (a, 6) that are quadratic non- residues having a ^ 1 and b ^ 1. 
Therefore 

#M,({a;,a;^a;n) = \{k -if -2 Q(A; - 1)) = \{k - l){k - 3). 

(Note that (1-1) is a quadratic residue, so the invertible non- residues 
of the form (a, 1) and (1,&) are disjoint.) A similar argument applies 
to the two remaining cases, where we rely on the fact that invertible 
elements of F^ x F^. and of F^. fall evenly into the six sextic residue 
classes. □ 
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The table in Proposition [33] gives the value of i^Mk{S) for the four 
subsets S* C /Xg that appear in Conjecture [301 It remains to construct a 
similar table for the values of jj^M]^\S). It turns out that these values 
can be expressed in terms of the number of points on a certain curve of 
genus four over various finite fields. We begin with a description of the 
curve that we need, after which we count points in order to compute 
the desired values. 

Proposition 34. Let ¥ he a perfect field of characteristic not equal 
to 2 or 3. For k G F* we define E^'^'> to he the elliptic curve 

and for 7, 5 G F* we define C^^'^^ to be a smooth projective model for 
the algebraic curve given hy the affine equation 

(a) The curve Cq^'^^ has genus four. 

(b) There are finite maps from Cq^'^^ to curves of the form E^'^^ given 
hy the following formulas: 

— ^ (x, z) ^ (-45x, 875/ - 45), 

6 ' ' " ^ ^ ^3 



^4 ' ,3 



^2 



(c) The maps in (b) are independent, hence they induce an isogeny 

(d) For any prime £ different from the characteristic of ¥, we have 
isomorphisms of Gal{¥/¥) -modules, 

= Hl,{Eff\Q,) X Hl,{Ef/''\Q,) 

xHl{E%'''\Q,)xHl{E%^''\Q,). 

Proof, (a) All of the C^'^^ curves are geometrically isomorphic, so it 
suffices to calculate the genus of C^'^\ which for convenience we de- 
note Ce. A simple calculation shows that the projective closure of 
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in is singular at (0, 0) and at the point at infinity, and that each 
of these singular points resolves to three points on the smooth model. 
(See Proposition [35] for details.) We let Ci be the elliptic curve 

Ci : z{l -z)= a;^ 

and we consider the natural degree 6 map 

The map ■?/' is ramified only at (0,0) and oo, the sets '?/'~^(0,0) and 
i\}~^{oo) each consist of three points, and each of these points has ram- 
ification index 2. Applying the Riemann-Hurwitz genus formula to i\} 
gives 

2(7(C6)-2 = 6(2(7(Ci)-2)+ 5^ (ep(V^)-l) = 6(2 - 2) +6(2 - 1) = 6. 

PeCi 

Hence g[C^ = 4. 

(b) It is an exercise to verify that the given maps are well-defined, but 
we briefly comment on their origin. The automorphism group of the 
curve Cp'^"^ is fairly large, since 

M3 X C Aut(C('^'')), [C,^]{x,z) = {Cx,^z). 

Taking quotients of Cq^'^^ by various subgroups of /X3 x fig gives maps 
to curves of lower genus, which in turn give the four maps described 
in (b). 

(c) From general principles, the maps in (b) induce isogenics 

4^'^^ for the given values of k. There are various ways to see that 
these isogenics are independent. For example, one can use the fact 
that the four E^'^^ are non-isogenous over C(7,5), treating 7 and 5 as 
indeterminates. Or, at least in characteristic 0, one can take 7 = 5 = 1, 
untwist to get four maps C^^'i) ^ defined over O, and use the ac- 
tion of Gal(Q/Q) on the maps to show that they are independent. (See 
AppendixlHl) Or, for a purely geometric proof, one can use intersection 
theory and the fact that the pairing 

(■,-):Map(C(^'^\E«)/^^'^-Z, 
(0, ip) = deg(0 + ip) - deg(f)- deg ^, 

is a positive definite quadratic form. (The E^^^ in the denominator is 
shorthand for the right action of the group of translations.) 

(d) It is a standard fact that H^^ of a curve and its Jacobian are isomor- 
phic. This gives the first isomorphism, and the second follows from (c) 
and the fact that an isogeny between abelian varieties induces an iso- 
morphism of their etale cohomologies. □ 
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Proposition 35. Let ^ be a prime ideal in Ok such that /Xg C Ok/ ^, 
i.e., Ni^/Q(^) = 1 (mod 6). Let ( E and C, G ^3, choose ele- 
ments 7, 5 G Ok satisfying (-^g = C o.'^^d {-^^ = ^, <^nd let C^''^'' be 
the smooth projective curve from Proposition given by the affine 
equation 

Then 

#MW(C,0 = ^(#^^r^ (^)-e(C,0 
where the error term e(C,0 given by the formula 



e(C,e) 



6 tfC = l 
^/C7^1 



+ 



3 ife = ^' 
^/cv^ 



+ 



3 ^/C' = e 

^/CVe 



Proof. Our choice of 7 and 6 imply tliat for any A G Ok, 

6 

= ^ (5~^A(1 — A) = non-zero cube (mod ^). 



A(l-A) 



c 



7 ^A = non-zero sixth power (mod 



We thus get a natural map 



'x, z) G a 



(7,<5) 



: X 7^ 0, 00 

(a;, 2;) 



M^(C,0, 



(31) 



We claim that the map 0311) is exactly 18-to-l. To see this, let A G 
Msi{(, Then A = 71;^ (mod ^) and A(l — A) = (mod for some 

■u, f G (C/</^)*, so A is the image of the point {u,v) G C^'^\Ok/^)- 
Further, for a given value of A, there are six choices for v and three 
choices for u. (Note that Ok/^ contains ^g.) Hence 



^(7:5) 



o 



K 



\ {x = or 00} 



It remains to count the number of /^-rational points with x = 
or 00 on a smooth model of C, 



(7,5) 



To ease notation, we let C = Cg'^'"^ and we let C be the curve 



.(7,5) 

C : 7(1 - iz') 



6x^ 



(32) 



The birational map 

C 



{x,z) 



[xz 
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is a bijection on the set of points 

C \ {x = or cx)} C" \ {x = or 00} U {2 = 0}, 

and the affine piece of C defined by equation (!32l) is smooth, so the 
points with x = on C become the points with x = or 2 = on C . 
(More precisely, we will see that the singular point (0, 0) G C is blown 
up to three points on C", while there are six smooth points of the form 
(0, 7~^/^) on both C and C .) The points on C with x = or 2; = 
are characterized by 

(0,2) eC'^ 2^=7"^ and (x, 0) e C" ^ x^ = 7r \ 

Thus there are points of the form (0, z) if and only if = 1, and 

there are points of the form (0,x) if and only if (-'^^jg = 1- Using the 

values = C and {^)^ = {^^H^), = CTS we find that 



# (0,;^)gC, 



(7,5) / ^ 
6 



# (x,0)gC, 




(7,5) / Ok 

It remains to count the points at infinity on C . Homogenizing the 
equation for C gives the curve 71/^ — 7^2;^ = Sx^y^. The unique (singu- 
lar) point at infinity is [x, y, z] = [1, 0, 0], so dehomogenizing by setting 
X = 1 gives the curve 

72/^ - = 5y\ 

We blow up the singular point (0, 0) by setting y = z^u. (This corre- 
sponds to blowing up twice. One can check that the other coordinate 
charts do not yield any additional points.) The resulting curve has 
affine equation 

-^z^u^ - 7^ = 5u^. 
This affine curve is smooth, and the points that map to the point at 
infinity on C are the points with 2; = and v? = —7^5^^. Using 
(^)3=(i)6(^)3'-^-«that 



# < points at infinity on C, 




(7,5) f Ok_ 

' V ^ 

This completes the proof of the proposition. □ 

The next step is to count the number of points on Cg'^''^^ defined 
over a finite field. This is done using the decomposition of Jq^'^"^ into a 
product of elliptic curves. 
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Proposition 36. With notation as in the statement of Proposition lS^ 
choose an element n G Ok satisfying A = hOk o-nd tt = 2 (mod 3). 
Further let e = (j)^- Then 

#4-^''^ (^) = ^^/Q + ^ + Tr^/Q(evr) + TTK/Q{e'C'e^) 



// ^ is an inert prime, say ^ = kOx with G Z satisfying k = 
2 (mod 3), and if we take 5=1, then the formula simplifies to 



k' + l + 8k z/C = l, 



P + l + 2k z/C^±l. 



Proof. To ease notation, let Fj^ = Ok /A, so Nk/q^ = i^^si- Further 
let be the (N^/q ^)*'^-power Frobenius map on F^^. Then the number 
of points in Ci^'^\¥si) is given by the trace formula [6l C.4.2], 

#C(^'')(F^) = Nk/q^+1 - Tr(F^ | Hl{Cij^lQi)). (33) 

We compute the trace using Proposition [311 which splits the represen- 
tation for Cg''''^^ into a product of representations on elliptic curves with 
zero j-invariant. Thus 



Tr(F^ I HiiC^ 



(7,5) 
6/Fjr 



= Tr (f^ I HUE%f\ Q,)) + Tr (f^ | H},{Ef;^''\ Q 

+ Tr(F^ I Hi{E''/;f\Q,))+TT[Fs, \ Hl{E^^;f\Qe)) ■ (34) 

We now apply [3, Chapter 18, Section 3, Theorem 4], which gives a 
formula for the trace in terms of residue symbols. Writing ^ = ttOk 
with TT = 2 (mod 3), we find that 



2652\-i /2V'5'\~' /2V5' 

vr — — — TT — — - — vr — — - — I vr 



vr — TT — TT — I TT 



*^ '6 V *^ / 6 V *^ / 6 *^ 
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-1 



C 4 TT 



2 \-^/ 3 



2 \-^/ 3 



(35) 



Noting that {^)^ = (-l)(Nif/Q^-i)/2^ combine §^ and dSS]) to 
obtain the desired resuh. 

In the case that ^ = kOx is an inert prime, we have (—1)2^*'' = 1 
since k is odd. Further, both 2 and k are primary, so cubic reciprocity 
gives = (1)3 = (1)3 ~ Further taking 5=1 imphes that ^ = 1, 
so the formula for ^Cq"''^\Ok/^) becomes 

k^ + l+ (Tr;^/Q(l) + Tr;^/Q(C') + Tr^/Q(C') + TTK/Q{0)k. 
Taking the six possible values ( & fi^ yields the stated formula. □ 

Proposition 37. Let k > 5 be a rational prime. The following table 
gives the values ofj^Ml\S,l) for various subsets S C fJ,Q, divided into 
cases according to whether k is split or inert in K = Q(-\/— 3), cf. 
Proposition 





k = 1 (mod 3) 


k = 2 (mod 3) 


(a) 




1(^-1)' 




(b) 




l{k-l){k-3) 


- 1) 


(c) 


#Mf'({c.,a;^cc;^ 6.^,1) 






(d) 


#Mr(M6,l) 


lie -2k + 4) 


lie + 2k -2) 



Proof. We begin with the case that k = 2 (mod 3), so = kOx is a 



prime ideal with Nk/q^ = k"^. We let uj = ^(1 + \/—3) be the usual 



sixth root of unity, and we choose some 7 G Ok satisfying 



u. 



Then for any < z < 5 we have 

18#mW(^M) = #C^'"''^M- 

from Proposition [351 with ( = ui^ and ^ = 1, 



12 if i = 
6 if i = 3 
otherwise 
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'k'^ + l + Sk if Z = 
k'^ + l-Ak if 2 = 3 
k"^ + 1 + 2k otherwise 



12 if z = 
6 if 2 = 3 
otherwise 



from Proposition [36] with ( = uj^ and ^ = 1, 

k'^ + 8k-ll ifi = 0, 
P - 4A; - 5 if i = 3, 
k"^ + 2k + 1 otherwise. 

It is now easy to compute ij^M]^^ (S*, 1) = Xl^es "f^^k^^ (C; 1) fo^ the four 
cases of the proposition. For example, 

#mW(^6,1) = ^{{e + 8k-n) + {e -Ak-5)+A{e + 2k + l)) 
= hk^ + 2k~2). 

o 

Next we consider the case that k = 1 (mod 3), so kOx = sphts. 
The definition of the residue symbol says that 



A 



A(l-Ar 
kOf 



A(l-A)\ A(l-A) 



-2v\ 



'K / 3 

SO using the Chinese remainder theorem 

a quantity such as m]1\c,,C) breaks up into a sum of products 
Hence for < i < 5 we have 

(For the second equality we've used the identity M^^ {CiCi = (C 0- 
We choose 7 and S to satisfy 

-^1 = u! and I — ; 



(36) 
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Then Proposition [35] gives us the formula 

ISMWK,^;^^) = i^Cf^'^\¥^) - e{uo\u'^), (37) 

where 

3 ii u = v (mod 3) 
otherwise 

3 if 2u = t> (mod 3) 
otherwise ' 

Further, Proposition [36] gives us the number of points on the curve, 

+ Tr^/Q(ecu5"+2.-) ^ Tr;,/Q(e^"+2.-)^ (gg) 

where e= (D3. 

Combining ( 136]) . ( [371) . ( l38l) . and ( [391) gives an explicit, albeit quite 
complicated, formula for m].^' (cj*, 1). In principle, this formula could be 
computed by hand, but we are content to evaluate it using PARI [26j, 
which yields the following values: 

{F + 4A; + 13 iii = 0, 
+ 1 if z = 3, 
fc^ — 2fc + 1 otherwise. 

(See Remark [38] for further information about this computation.) It is 
now a simple matter to compute the value of #m|^' (5, 1) for the four 
cases. For example, 

#mW (^g, 1) = 4 ((fc' + 4fc + 13) + (fc2 - 8A; + 7) + 4(fc2 - 2fc + 1)) 
18 

This completes the proof of Proposition [37] □ 

Remark 38. We used PARI [26j to compute m]!i\uj^, 1) by evaluating 
formulas (l36l) . (137]) . (138]) . and (1391) . where we treated fc, tt, and e as 
indeterminates, and we formally set vf = fc/vr and e = 1/e. (The 
PARI code used for the computation is given in Appendix [C]) The 
value of m]1\uj'^, 1) turns out to be a quadratic polynomial in k that is 
independent of k mod 4. We do not have an a priori explanation for 
why this should be the case. In order to illustrate the dehcacy of the 



e(u;",u; 



2v\ 



6 if M = (mod 6) 
otherwise 



+ 
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argument, we suppose for a moment that the isogeny decomposition of 
the Jacobian of C^"''^^ in Proposition [31] looks Hke 

(All that we have done is change the second elliptic factor from _E(^'>'^<5*') 
to E^^^'^\) This would have the effect in formula (1391) of changing the 
second trace term from Trii:/Q(e^u;'^""'"''^7f) to Trx/Q(e^ci;^^7r). But with 



this small modification, there is less cancelation in the computation 
of m|.^^ {uj\ 1), so for example ij^M]^^ ({u;, u;^}, l) would equal 

Thus ij^Mj^'^ [{uj, 00^}, l) would depend on both and on the factor- 
ization of k in Ok- 




Remark 39. Many of the cases of Proposition [371 can be obtained some- 
what more easily by working on elliptic curves z{l — z) = Sx^ or genus 
two curves 72:^(1 — 7-^^) = Sx^. However, some cases require the curves 
7z^(l — 7-2®) = Sx^ of genus four, so for unity of exposition and to save 
space, we have derived all cases using these latter curves. 



Combining Conjecture [30] with the computations in Propositions 
and [37] yields precise formulas for the conjectural density of Type 1 
primes on y 



•"^ — x^ + k when k is prime. 



Conjecture 40. Let k > 5 be a rational prime. Then 



lim 



- + R{k), 



where R{k) depends on k (mod 36) and is given by the following table: 





k mod 36 


R{k) 






k mod 36 


R{k) 


(a,c) 


1, 19 


2 




(a,c) 


17, 35 


2 


3(A;-3) 




3ik-l) 


(b) 


13, 25 







(b) 


5, 29 





(d) 


7, 31 


2k 




(d) 


11, 23 


2k 


3{k-2y 




3{k^-2) 



In particular, R{k) = 0{l/k). 
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We do not have an intrinsic explanation for why R{k) is the same in 
cases (a) and (c), nor do we know why R{k) = in case (b). 

7. Amicable pairs for elliptic curves — Experiments 

In this section we present the results of experiments that test the 
reasonableness of our conjectures. We begin with Conjecture [3l which 
deals with the case of CM curves having nonzero j-invariant. 

We computed the number Qe{X) of amicable pairs up to X for 
elliptic curves with CM by the imaginary quadratic order of discrimi- 
nant —D and conductor /. Theorem [13] says that it suffices to consider 
D = 3 (mod 4). Further, the assumption that E is defined over Q 
means that O has class number one, so the classification of imaginary 
quadratic fields of class number one combined with an elementary for- 
mula for the class number of an order |18l Exercise 4.12] imply that 
the only possibilities for D are D G {3,7,11,19,43,67,163}, and the 
possible values of / are given by / G {1,2,3} if D = 3, / G {1,2} if 
D = 7, and / = 1 in all other cases. See [IHl A §3] for a Weierstrass 
equations for each CM type. 

We ignore for the moment the case (-D, /) = (3,1). As noted in 
the proof of Theorem [T3l the curves with {D, f) equal to (3, 2), (7, 1), 
and (7, 2) have nontrivial 2-torsion, so neither they nor any of their 
(necessarily quadratic) twists have amicable pairs. The curve with 
{D, /) = (3, 3) listed in [TH A §3] has nontrivial 3-torsion, but it has 
quadratic twists with trivial torsion, so is a candidate to have amicable 
pairs. Table [2] lists the number Qe{X) of amicable pairs up to the given 
bound and the ratio of Qe{X) to the number Me^X) of primes p such 
that ^Ep(¥p) is prime. For this table we used the following Weierstrass 
equations^ 



(A/) = 


(3,3) 


y' 


= a;3 


- 120x + 506, 


(A/) = 


(11,1) 


y^ + y 


= x' 


-x^ -7x + 10, 


(A/) = 


(19,1) 


y^ + y 


= 


- 38x + 90, 


(A/) = 


(43, 1) 


y^ + y 


= 


- 860x + 9707, 


(A/) = 


(67, 1) 


y^ + y 


= x^ 


- 7370a; + 243528, 


(A/) = 


(163,1) 


y^ + y 


= x' 


- 2174420X + 1234136692 



The results in Table|5]are consistant with ConjecturedHl which predicts 
that the ratio QE{X)/AfE{X) should approach i. 



Calculations on quadratic twists of the listed curves yielded virtually identical 
results. 
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(A/) 


(3,3) 


(11,1) 


(19,1) 


(43,1) 


(67,1) 


(163,1) 


Qi?(l05) 


124 


48 


103 


205 


245 


395 


Qe(105)/Ar^(105) 


0.251 


0.238 


0.248 


0.260 


0.238 


0.246 


Qi^(106) 


804 


303 


709 


1330 


1671 


2709 


Q^(106)/Ars(106) 


0.250 


0.247 


0.253 


0.255 


0.245 


0.247 


Qi.(10^) 


5581 


2267 


5026 


9353 


12190 


19691 


Qs(107)/Ars(107) 


0.249 


0.251 


0.250 


0.251 


0.250 


0.252 



Table 2. Qe(X) and Qe{X)IMe{X) for elliptic curves 
with CM by Q(v^^) 



k 


2 


3 


5 


6 


7 


10 


X = 10^ 


0.251 


0.122 


0.081 


0.134 


0.139 


0.125 


X = 10^ 


0.250 


0.139 


0.083 


0.142 


0.133 


0.107 


X= 10^ 


0.249 


0.139 


0.082 


0.139 


0.129 


0.107 



Table 3. QE{X)/AfE{X) for elliptic curves y'^ = + k 



We next considered the curves y"^ = + k with j{E) = 0. Table [3l 
which is included for historical reasons, was our first intimation that the 
limiting value of Qe{X)/N'e{X) behaves differently for different values 
of k, with no obvious pattern for 2 < k < 10. (Note that we do not 
list values of k that are squares or cubes, since in those cases -^(0)1018 
is nontrivial, so there are no amicable pairs.) 

We recall the notation A/"]^^ for the set of Type 1 primes for the 
curve y'^ = x^ + k; see Section[6]for the precise definition. Conjecture [271 
predicts that Qfc(X) ~ J^ft\x), and in the case that k is prime. 
Conjecture HQ] says that 

a4''W~ Q + i?w)A4.(x), 

where R{k) is given by an explicit formula that depends on k mod- 
ulo 36. We tested these two conjectures by computing Qk{X), A/"]^^(X), 
and Mk{X) for X = 10*. The results are listed in Table HI Col- 
umn 5 provides convincing evidence for Conjecture [271 cind the final 
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two columns show that Conjecture HO is in good agreement with exper- 
iment in all eight cases. (The notation {x.n) after each value of k indi- 
cates the case x = (a),. . . ,(d) and the congruence class k = n (mod 3) 
from Conjecture HOI) 





Density of Type 1 primes 

Ml\x)IMu{X) 


k 


Sfe(X) 


ArW(X) 

k ^ ' 


A^fc(^) 




experiment 


coniecture 








-L 1 O 1 UO 






1 — n QQQQ 


7 (A A\ 


/I QQQ9 


o 1 oZO 


lDc5 1 "rtO 


n 9^1 
u. zo± 


u. ozuo 


— n 1^900 

— U-OZUU 


11 (d 2) 


'^'^791 




-Luyuuz 


u. zoo 




— n "^Qf^n 
YYg — u.oyou 


{u ^^ 






-LU ( OOO 


252 




1 _ fi QQQQ 

2^ — u.oooo 




OZUUo 


DOoiU 




n 9f^i 
u. zo± 


U.O M i 


— n 'i7'^n 

g — U-O 1 ou 


19 (c.l) 


31729 




1 RSI Q(S 


0.252 


0.3750 


3 _ n Qyc^fi 

g^ U .iJ 1 tJU 


23 (d.2) 


30480 


61210 


168512 


0.249 


0.3632 


# = 0-3624 


29 (b.2) 


28085 


56286 


168642 


0.249 


0.3338 


i = 0.3333 


31 (d.l) 


30301 


60349 


168344 


0.251 


0.3585 


ir = 0.3579 


37 (a.l) 


29728 


59430 


168471 


0.250 


0.3528 


^ = 0.3529 


41 (b.2) 


28050 


56381 


168567 


0.249 


0.3345 


\ = 0.3333 


43 (d.l) 


29619 


58807 


168410 


0.252 


0.3492 


W = 0-3504 


47 (d.2) 


29220 


58400 


168365 


0.250 


0.3469 


^ = 0.3475 


53 (a.2) 


29278 


58257 


168353 


0.252 


0.3460 


^ = 0.3462 


59 (d.2) 


29378 


58422 


168783 


0.252 


0.3461 


3479 ~ ^.0440 


61 (b.l) 


28027 


55816 


168197 


0.251 


0.3318 


\ = 0.3333 


67 (d.l) 


29242 


57944 


168239 


0.253 


0.3444 


= 0.3439 


71 (c.2) 


28789 


57661 


168508 


0.249 


0.3422 


11 = 0.3429 


73 (a.l) 


28975 


57828 


168614 


0.251 


0.3430 


i| = 0.3429 


79 (d.l) 


29127 


57937 


168690 


0.252 


0.3435 


^0^9 — n S422 
5929 — 


83 (d.2) 


29032 


57871 


168435 


0.251 


0.3436 


2351 _ n .14^1 4^ 
6887 - U.04i4 


89 (a.2) 


28909 


57634 


168737 


0.251 


0.3416 


II = 0.3409 


97 (b.l) 


28014 


55880 


168457 


0.251 


0.3317 


i = 0.3333 



Table 4. Density of Amicable and Type 1 primes with 
p < X = 10^ for the curve = + fc, prime k. 



We also checked Conjecture EOl experimentally for composite values 
of k. The results are listed in Table where the conjectural limiting 
ratio is obtained by explicitly counting the size of the sets M.}^ and ■ 
The top eight k entries in this table are products of two primes covering 
the usual eight cases; the final four entries include two values that are 
not square- free (175 = 5 ■ 7^ and 245 = 5-7^) and two values that are 
products of three distinct primes (385 = 5 ■ 7 ■ 11 and 455 = 5 ■ 7 ■ 13). 
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In order to further test the vahdity of Conjecture [30], we recomputed 
the final entry in the table with X = 10^ and obtained 

A45!5(10')/Ar455(10^) = 0.3380. 
This is in excellent agreement with the theoretical value of -^1^ = 

° 13915 

0.3377. 





Density of Type 1 primes 
Afl^\x)/jVk{X) 


k 


Qk(X) 






Q/AfM 


experiment 


conjecture 


35 (d.2) 


31423 


63169 


168666 


0.248 


0.3745 


= 0.3739 


55 (d.l) 


29645 


58718 


168870 


0.253 


0.3477 


Mf = 0-3467 


77 (b.2) 


28170 


56251 


168921 


0.251 


0.3330 


I = 0.3333 


85 (b.l) 


28187 


56142 


168767 


0.251 


0.3327 


i = 0.3333 


323 (c.2) 


28396 


56609 


168585 


0.251 


0.3358 


= 0.3359 


629 (a.2) 


28210 


56269 


168042 


0.251 


0.3349 


3267 _ r, QQAr 

9766 ~ U..3d40 


703 (c.l) 


28558 


56754 


168817 


0.252 


0.3362 


3278 ~ 0.3347 


901 (a.l) 


28341 


56384 


168411 


0.252 


0.3348 


3738 _ n no/i-i 
11189 U.OO-il 


175 (d.l) 


31543 


63177 


168840 


0.250 


0.3742 


4rF = 0.3739 

115 


245 (b.2) 


29722 


58848 


175934 


0.253 


0.3345 


i = 0.3333 


385 (b.l) 


28070 


56158 


168393 


0.250 


0.3335 


i = 0.3333 


455 (d.2) 


28346 


56627 


168342 


0.250 


0.3364 


4699 _ Q 3377 
13915 ~ ' 



Table 5. Density of Amicable and Type 1 primes with 
p < X = 10^ for the curve = + k, composite k. 



Finally, we consider Conjecture [9], which deals with non-CM curves. 
This conjecture is much harder to check numerically because the func- 
tion y/X/{\ogXy grows quite slowly. We performed an extended 
search for amicable pairs on the elliptic curve 

E : + y = + x'^ (40) 

of conductor 43 that we studied in Example [U Appendix [D] lists all 
normalized amicable pairs {p,q) with p < 10^^. Conjectures El^a) and [HI 
says that Qe{X), the number of amicable pairs up to X, should grow 
like a multiple of ^/X/{\ogXy. Table [6] tests this conjecture by com- 
puting the ratios 

Qe{X) , hgQEiX) 

—F= and — ; — — 

v^/2(logX)2 logX 

for various values of X. The computations for Table [6l i.e., the compu- 
tation of Q£;(10^^), took approximately five days running in parallel on 
a cluster of ten machines. The third column of Table [6] provides some 
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small support for the conjecture that Q{X) grows like a multiple of 
V^/(logX)^. On the other hand, although the fourth column of the 
table suggests that Q{X) grows like for some 5 > 0, it is far from 
clear that 6 is as large as |. We suspect the problem is that we are only 
able to compute Q{X) up to X = 10^^, and although 10^^ is a moder- 
ately large number in terms of computation time, it is comparatively 
small compared to the likely error terms in any putative asymptotic 
formula for Q{X). 



X 


Q{X) 




log S(X) 
logX 


10^ 


2 


0.382 


0.050 


10^ 


4 


0.329 


0.086 


108 


5 


0.170 


0.087 


109 


10 


0.136 


0.111 


lOlo 


21 


0.111 


0.132 


2 • IQio 


32 


0.127 


0.146 


4 • IQio 


37 


0.110 


0.148 


6 • IQio 


44 


0.111 


0.152 


8 • IQio 


53 


0.118 


0.158 


lOii 


55 


0.112 


0.158 



Table 6. Counting amicable pairs for y"^ + y = + x'^ 



Finally, we searched for normalized aliquot triples {p, q, r) on the 
curve (HOll . We found no examples with p < 10^. To compare this to 
Conjecture [HI we consider the counting functions 

[^] , , J, f/ \ Pi distinct primes with Pi < X and 

(X) = # . . . : ^^^^^^^^^ ^ ^^^^ for 1 < . < £ - 1 

Note that P^'(X) is not counting ahquot cycles, because we do not re- 
quire that ^Ep^(¥p^) equal pi. A natural generalization of the Koblitz- 
Zywina conjecture (Conjecture [6]) is that 

We computed V'e\X) and V%\X) for the curve ( Hoi) and used it to 
estimate the values of the constants C^3in and C^^m- The results are 
listed in Table [3 
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X 




V%\X) 








l& 


485 


21 


0.643 


0.320 


0.033 




3099 


116 


0.592 


0.306 


0.058 


10'^ 


22328 


741 


0.580 


0.310 


0.117 


10« 


168611 


4995 


0.572 


0.312 


0.250 



Table 7. An aliquot triple estimate for + y = + 



Our heuristic argument from Section [3] suggests that if (p, g, r) is a 
triple counted in (X), then the probability that ^Er(¥r) equals p 
is 0(r~^/^), and more precisely, we expect it to be between 
and ir~^/^. Thus the number of aliquot triples less than X should 
roughly be somewhere between 

IV^^\X)X~'/' and \vf{X)X-'/\ (41) 

The last column of Table [7] gives our heuristic upper bound for the 
number of aliquot triples up to X. It is thus not surprising that we 
found no examples up to 10^. Using the estimate C^^/q ~ 0.312, setting 
jj^^xjs, and taking the upper bound in (jH]), it would 
require at least X ^ 5 ■ 10^ in order to expect find an aliquot triple. 
And even for X = 10^^ we wouldn't expect more than two or three. 

8. Motivation and generalizations 

Remark 41. Elliptic amicable pairs and aliquot cycles appeared in a 
natural fashion when the authors were generalizing to elliptic divisi- 
bility sequences Smyth's results ^3J on index divisibility of Lucas se- 
quences. Let {Wn)n>i be a normalized nonsingular nonperiodic elliptic 
divisibility sequence associated to an elliptic curve E/Q, and consider 
the set 

S = {n>l:n\ Wn}. 

An element of S is called basic if n/p ^ S for all primesp dividing n. 
It turns out that basic elements can be created using aliquot cycles, 
a phenomenon that does not occur in Smyth's work. More precisely, 
given any aliquot cycle {pi, ■ ■ ■ ,pe) for E/Q of length £ > 2, the prod- 
uct PiP2---Pe is a basic element for (Wn), and more generally, any 
product of such products is basic. See [H] for details. 

Remark 42. As we have defined them, aliquot cycles for elliptic curves 
differ in a significant way from classical aliquot cycles associated to the 
sum of divisors function. In the classical case, every integer n leads to 
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a possibly non-repeating aliquot sequence (n,a{n),a'^{n),a^{n), . . 
and it is an aliquot cycle if some iterate (T^{n) eventually returns to n. 
(A major open problem is whether there are starting values for which 
the sequence is unbounded.) But for elliptic curves, if we arrive at 
a prime p such that jj^Ep{¥p) is not prime, then the sequence cannot 
be continued. We propose here two alternative definitions of elliptic 
aliquot sequences that more closely resemble the classical definition. 
We leave the investigation of these generahzed sequences to a future 
paper. 

Definition. Let -E/Q be an elliptic curve, let L{E/Q, s) = ^„>;^ o,n/n^ 
be the L-series of E, and define a function 

Fe-.N^N, FE{n) = n + l-an. 

A Type L aliquot sequence for E/Q is a sequence obtained by starting 
at some n G N and repeatedly applying the map Fe- A Type L aliquot 
cycle is a Type L aliquot sequence that repeats at its starting value. 

Definition. Let E/Q be an elliptic curve, let S^/Z be the open subset 
of the Neron model for E/Q consisting of the connected components 
of each fiber, and define a function 

Gij : N — ^ N, GE{n) = #^°(Z/nZ). 

A Type N aliquot sequence for i?/Q is a sequence obtained by starting 
at some n eN and repeatedly applying the map Ge- A Type N aliquot 
cycle is a Type aliquot sequence that repeats at its starting value. 

Remark 43. There is a natural way to generalize the notion of elliptic 
amicable pairs and aliquot cycles to elliptic curves defined over number 
fields. Thus let F/Q be a number field and E / F an elliptic curve. We 
will say that a sequence of distinct degree one prime ideals pi, • • ■ , p£ 
is an aliquot cycle of length ^ for E/F if E has good reduction at 
every pj and 

#Ep,(FpJ = Nx/q(P2), #^p.(FpJ = N^/q(P3), . . . 

#Ep,_,(Fp,.J = N^/q(p,), #^p,(FpJ = Nk/q(Pi). 

Many of the methods and results in this paper carry over in a straight- 
forward manner to the number field case. For example, the following 
analogue of Theorem [13] holds. 

Theorem 44. Let F/Q he a number field, and let E/F he an elliptic 
curve with complex multiplication by an order in the quadratic imag- 
inary field K. Suppose that p and q are degree one primes of F at 
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which E has good reduction, that ISlp/Qp > 5, and that 

#4(Fp)=.N^/Qq. 

Assume further that j{E) ^ 0. Then 

#E,{¥,) = N^/Q p or #^q(Fq) = 2 N^/q q + 2 - N^/q p. 

It would be interesting to see to what extent tlie otiier results in this 
paper are valid over number fields, including especially the analysis of 
amicable pairs on curves with j{E) — 0. 
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Appendix A. Curves with j = have no aliquot triples 

In this section we use Corollary [24] and a detailed case-by-case anal- 
ysis to show that an elhptic curve with j = has no normalized aliquot 
triples (p, q, r) with p > 7. The details are sufficiently intricate that it 
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seems likely a different argument would be needed to prove that there 
are no aliqout cycles of length greater than three. 

Proposition 45. Let E/Q be an elliptic curve with j{E) = 0. Then E 
has no normalized aliquot triples [p, q, r) with p > 7. 

Proof. We use Corollary [24l which says that if p and q = j^Ep(¥p) are 
prime, then r = ^Eg(¥g) takes one of six possible values. One of these 
six possible values is p, which is not allowed since we are assuming 
that p, q, r are distinct. Hence r has one of the following forms, 

r = 2g + 2 - p, (Case 1) 



±iq + l-p)± 3Ap,, 



(Case 2) 



where A^-^y satisfies 



2 4:xy — [x + y 



A 



(Of course. Case 2 is really four cases, depending on the choice of signs.) 

For the moment letting s = ^Er{¥r), we can apply the same rea- 
soning to (g, r, s) to deduce that 

s = 2r + 2 - g, (Case A) 



±(r + 1 - g) ± 3Ag^r 



2 

To ease notation, we let 

F{x,y) = 



(Case B) 



±{y + l-x)±3A,^y 



2 

Then the two cases for r followed by the two cases for s give four 
possibilities for s in terms of p and q: 

s = 3q + Q-2p, (Case lA) 

s = 2F{p, q) + 2-q, (Case 2A) 

s = F{q,2q + 2-p), (Case IB) 

s = F{q,F{p,q)). (Case 2B) 

(Of course, each case is really several cases depending on the choice of 
signs for each occurrence of F.) 

The assumption that [p, q, r) is an aliquot triple is equivalent to 
saying that s = p. Suppose first we are in Case lA. Then s = p is 
equivalent to 

3g + 6 — 2p = p. 
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SO p = q + 2. This contradicts our assumption that the triple is normal- 
ized, i.e., that p is the smallest element of the triple. Hence Case lA 
is not possible. 

Next we consider Case 2A. Then the assumption s = p implies that 
2F{p, q) = p + q — 2. Using the definition of F, this can be written as 

±{q + 1 - p) ±3Ap,g = p + q - 2, 

which (using the definition of A) implies that 

((p + g - 2) ± (g + 1 - p))' = 9Al^ = 3(4pg - (p + g - l)^). (42) 

This gives two subcases, which we denote by 2A"'" and 2A~ according 
to the choice of sign. A little bit of algebra yields 

28p2 - 24pg + 12g2 - 72p - 24g + 48 = 0, (Case 2A+) 

12p2 - 24pg + 28g2 - 24p - 40g + 16 = 0. (Case 2A-) 

Both of the functions on the left-hand sides have leading quadratic 
forms that are positive definite, so there are only finitely many integral 
solutions (p, q). A more careful analysis shows that the first is positive 
for p > 5 and the second is positive for p > 7. 

Next comes Case IB, where the assumption that s — p leads to the 
formula 

F{q,2q + 2-p) = p. 

Writing this out in terms of Aq^2q+2-p, moving all the other terms to the 
other side, squaring, and simplifying, we again get two cases depending 
a choice of sign. Thus 

12p2 - 12pg + 4g2 - 24p +12 = 0, (Case 1B+) 

4p2 - 4pg + 4g2 + 12 = 0. (Case IB") 

The quadratic function for Case 1B+ is positive for p > 7 and the 
quadratic function for Case 1B~ is positive for p > 0. 

Finally we turn to Case 2B, which is somewhat more complicated 
because it is given by the formula 

F{q,F{p,q)) = p, 

which involves two iterations of the function F. The signs on the A^^y 
terms arc irrelevant since we square them, but the other signs in the 
definition of F do affect the eventual equation. After a bunch of alge- 
bra, we find that the p and q values for an amicable triple coming from 
Case 2B must satisfy one of the following equations. 

4p^ + 2p^q + 3p^q^ - pq^ + q"^ - 6p^ - 15p^q 

- Ihpq^ + 3p^ + 3pg + = (Case 2B++) 
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9p2g2 - 9pq^ + 9q^ + 9p^q - 27pq^ + 3p^ - 21pq 

-3g2-6j9 + 6g + 4 = (Case 2B+-) 
3p2g2 - 3pq^ + q^ + 9p^q - 9pq'^ + 9p^ - 9pq + 3q^ = (Case 26^+) 
V - 18j9^g + 33p2g2 - 27pq^ + 9q^ - lOp^ + 33p^q 

- 2lpq^ + 2 V - 21j9g - 3g^ - 10j9 + 6g + 4 = (Case 2B"") 

All of these quartic functions are positive if < p < g with p sufficiently 
large. More precisely, it suffices to take p > 3 for Cases 26++ and 2B''"", 
p > 4 for Case 2B~+ and p > 2 for Case 2B . 

This completes the proof that E has no aliquot triples. □ 



Appendix B. Proof of independence of the maps in 

Proposition 



We give a proof for char(F) = 0. Independence of maps is geometric, 
so it suffices to prove independence for 7 = 5 = 1. Let ■i/'i, . . . , '?/'4 be 
the four maps in Proposition [311(b), so 

We compose these maps with untwisting maps so we get 

four maps 

(x,y)^{ \^ x,\y ] 

01 : Ci'''^ ^ ^(16) ^ ^ A ^(1)^ 

03:^(^-1) ^ i^(i) i^'y)-(-'y\ Eii)^ 
0^ : ^(-1) (-'y)^(-'^y), 

The maps ■j/'i , . . . , '?/^4 are defined over Q, but the maps 0i, . . . , 04 are 
only defined over Q, not Q. We consider the action of Gal(Q/Q) on 
these maps. To do this, we choose elements cr, r G Gal(Q/Q) satisfying 

fr(v^) = pv^, a{{} = i, 

r(^) = ^, r(z) = -z. 

Here p = |(— 1 + a/— 3) is a fixed primitive cube root of unity. We 
also note that /X3 acts on E^^^ via [p\{x,y) = {px,y). Looking at the 
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explicit formulas for 0i, . . . , 04, we find that 

0^=[p2]o0i, 0-=[p]o02, 03=' 



-1] O 04 



Now suppose that we have a relation 

[rii] O 01 + [n2] O 02 + [ns] O 03 + [m] O 04 = 0. (43) 

Applying the transformation r to ( l43i) has the effect of replacing 04 
by [—1] o 04, so subtracting the two equations yields [2n4] o 04 = 0. 
Since the map 04 is a finite map, it follows that 77,4 = 0. 

Applying a and to fj43|) . we end up with three equations 

hi] O 01 + [^2] O 02 + [^3] O 03 = 0, (44) 

[ni] o O 01 + [^2] o [p] o 02 + [7^3] O 03 = 0, (45) 

[ni] o [p] o 01 + [^2] o [p^] o 02 + [n'i] o 03 = 0. (46) 

Adding (j44]), (j45]), and(|46l) and using l + p + p^ = g gives [8^3] o 03 = 0, 
which implies that 77,3 = 0. Similarly, adding (jHj) to [p] times ( H5l) 
to [p^] times (H6l) gives [3ni] o 0i = 0, so ni = 0. Finally, since rii = 
ris = 0, the equation (l44l) gives 722 = 0. This completes the proof 
that 01, . . . , 04 are independent. 

Appendix C. Evaluating the formula for m|.^^ when 

k = 1 (mod 3) 

In this appendix we give the PARI f26] script that we used to com- 
pute Mj^\s, 1) via the formulas ((361), ([37]), (I3HD, and ([39]). In these for- 
mulas we treat k, tt, and e as indeterminates and formally set 7f = fc/vr 
and e = 1/e. The value of Mj^\S, 1) turns out to be independent of 
k mod 4 and is always a quadratic polynomial in k. The output from 
the routine TestFormulaForSplitPrimes is given in Table [8|, which 
was calculated using the following PARI program. 

/* The function TestFormulaForSplitPrimes computes M_k"{[l]} 

for k = 1 (mod 3), using the formula as a sum of products of #C - e terms. 
The following global variables must be left as indeterminates: 
k, pi, e2 

Here e2 represents the cubic residue of 2 modulo pi. 
The conjugates of these are given by 
pibar = k/pi and e2bar = l/e2. 

Further, w is assigned the value quadgen(-3) , and wbar = conj (w) = 1/w 

*/ 

{ 

TestFormulaForSplitPrimes () = 
local (m, ISets , IndexSet) ; 
print("\\begin{align*}") ; 
forstep(kmod4 = 1, 3, 2, 
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k = l (mod 4) 








= (l/18)(fc' 


+ 4A; + 13) 




= (l/18)(fc' 


- 2A; + 1) 




= (l/18)(fc' 


-2k + l) 


#MW(a;3^1) 


= (l/18)(fc' 


-8k + 7) 


#A^l(a.^l) 


= (l/18)(fe' 


-2/c + l) 




= (l/18)(fc' 


- 2fc + 1) 



/c = 3 (mod 4) 

(l^°, 1) = (1/18) (fe2 + 4fe + 13) 
#MW(^\l) = (l/18)(fe2-2A: + l) 
#MW(a;2^1) = (l/18)(A:^-2fc + l) 
#MW(a;^l) = (l/18)(fc2-8fc + 7) 
#MW(a;^l) = (l/18)(A:2-2A: + l) 
#MW(a;5, 1) = (1/18)(A;2 - 2A; + 1) 

#mW(V,c.5},i) = (l/9)(fc2_2fc + l) 
#mW(VScu3,c.5},i) = (l/6)(fc2 _4A: + 3) 
#MW({a;^a;^a;^a;5},l) = (2/9)(A;2 - 2fc + 1) 
#MW({a;^a;^a;^a;^a;^a;^}, l) = (1/3)(A;2 - 2A; + 4) 

Table 8. Results of computing of M^^\s, 1) using PARI 



printC kWequiv ",kmod4," \\pmod{4} WW"); 
for (i = 0, 5, 

m = MMSumC [i] ,kmod4) ; 

printK" W#M_k-{[l]}(Wo-",i,",l) &= "); 

prlntK" (" .content (content (m) ),")(" ,m/content (content (m)) ,")"); 
if (i < 5 I I kmod4 == 1, printC WW"), print); 

); 

if (kmod4 == 1, print (" WW")); 

); 

print("Wend{align*}") ; 
print("Wbegin{align*}") ; 

ISets = [[1,5], [1,3,5], [1,2,4,5], [0,1,2,3,4,5]]; 
for (j = 1, #ISets, 
IndexSet = ISets [j]; 
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m = MMSumdndexSet.l); 

printK" \\#M_k~{[l]}\\bigl(\\{"); 

for (i = 1, #IndexSet, 

print 1 ("\\o"" , IndexSet [i] ) ; 

if (1 < #IndexSet, printK",")); 

); 

printl("\\}, l\\bigr)\n k= (", content (content (m) ),")") ; 
printK" (" .m/content (content (m) ),")"); 
if (j < #ISets, print(" WW"), print); 

); 

print("\\end{align*}") ; 

} 

w = quadgen(-3) ; 
wbar = conj (w) ; 
pibar = k/pi ; 
e2bar = l/e2; 

{ 

CC(u,v,kmod4) = 

if (kmod4 == 0, error("k mod 4 must be 1 or 3")); 
k + 1 

+ w" (2*v) *pibar + wbar" (2*v) *pi 

+ e2"2*w"(3*u+4*v)*pibar + e2bar"2*wbar~ (3*u+4*v) *pi 
+ e2*w"(5*u+2*v)*pibar + e2bar*wbar" (5*u+2*v) *pi 
+ (-!)"( (kmod4-l) /2) *e2*w" (u+2*v) *pibar 
+ (-!)"( (kmod4-l) /2) *e2bar*wbar" (u+2*v) *pi ; 

} 
{ 

ee(u,v) = 

if (u 7. 6 == 0, 6, 0) 

+ if((u-v) ■/. 3 == 0, 3, 0) + if((2*u-v) "/, 3 == 0, 3, 0); 

} 

HM(u,v,kmod4) = (1/18) * (CC(u, v,kmod4) - ee(u,v)); 
{ 

MMSujii(IndexSet,kmod4) = 
local (i , s) ; 
s = 0; 

for (j = 1, #IndexSet, 
i = IndexSet [j] ; 
for (u = 0, 5, 
for (v = 0, 2, 

s = s + MM(u,v,kmod4)*MM(u-i,v,kmod4) ; 

); 
); 

); 

return(s) ; 

} 

Appendix D. Amicable pairs for y"^ + y = + up to 10^^ 

We used PARI-GP [26] to compute all normalized amicable pairs 
(p, q) on the curve y"^ + y = x"^ + x"^ with p < 10^^. The list is given in 
Table 1 
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(853,883) 
(1147339,1148359) 
(82459561,82471789) 
(253185307,253194619) 
(794563993,794571803) 
(2185447367,2185504261) 
(4101180511,4101190039) 
(5293671709,5293749623) 
(7074693823,7074704971) 
(9395537549,9395559011) 
(9849225103,9849306373) 
(12657210407,12657303353) 
(13789895011,13790023199) 
(14976551207,14976590371) 
(15679549877,15679688491) 
(17725049203,17725142719) 
(31615097957,31615194739) 
(33963999907,33964128017) 
(39287748091,39287808559) 
(46438194193,46438453213) 
(51881025571,51881167549) 
(55823622193,55823919169) 
(62765305697,62765625749) 
(66252308051,66252349753) 
(69449506103,69449741239) 
(77264683829,77264993327) 
(79067605783,79067881429) 
(94248260597,94248586591) 



(77761,77999) 
(1447429,1447561) 
(109165543,109180121) 
(320064601,320079131) 
(797046407,797057473) 
(2382994403,2383029443) 
(4686466159,4686510971) 
(6677602471,6677694539) 
(7806306133,7806380963) 
(9771430993,9771433303) 
(10574564857,10574619851) 
(13003880317,13003900901) 
(14436076927,14436180091) 
(15597047659,15597075937) 
(16322301811,16322366867) 
(17841395323,17841406601) 
(33266376239,33266419807) 
(34525477799,34525684639) 
(40136806357,40137038941) 
(51838270219,51838493561) 
(52011956957,52012184953) 
(57920520199,57920640709) 
(62995853671,62996152237) 
(67177409329,67177631771) 
(75002612911,75002660263) 
(77635421531,77635670141) 
(81263083703,81263204563) 



(94248260597,94248586591) | 

Table 9. Amicable pairs for + y = + up to 10^^ 



Mathematics Department, Box 1917, Brown University, Providence, 
RI 02912 USA 

E-mail address: jhsOmath . brown . edu 

Department of Mathematics, Simon Fraser University, 8888 Uni- 
versity Drive, Burnaby, BC, Canada V5A 1S6, and Pacific Institute 
FOR THE Mathematical Sciences, 200 1933 West Mall, Vancouver, BC, 
Canada V6T 1Z2 

E-mail address: kestange@pims.math.ca 



